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ABSTRACT 

A  Correlative  Comparison  of  Geomagnetic  Storms  and 
Aurorc!  Substorms  Using  Geomagnetic  Indices 

by 

William  B.  Cade  III,  Master  of  Science 
Utah  State  University,  1993 

Major  Professor:  Dr.  Jan  J.  Sojka 
Department:  Physics 

Geomagnetic  storms  and  auroral  substorms  are  usually  defined  in  terms  of 
their  respective  current  systems  and  measured  by  geomagnetic  indices.  While  there 
should  be  a  relationship  between  the  two  current  systems,  finding  one  has  eluded 
most  researchers.  However,  when  a  time-weighted  accumulation  technique  is  applied 
to  the  auroral  electrojet  index,  a  clear  relationship  emerges  which  shows  that  the  two 
current  systems  are  indeed  dominated  by  the  same  directly  driven  process,  differing 
only  in  how  they  respond  to  the  same  input.  Furthermore,  the  possibility  is 
demonstrated  of  modeling  and  even  forecasting  the  ring  current  firom  time-weighted 
electrojet  data. 


(118  pages) 


CHAPTER  1 


INTRODUCTION 

One  of  the  first  comprehensive  theories  of  geomagnetic  storm  development 
was  put  forth  by  Chapman  and  Bartels  [1940].  They  assumed  that  most  of  the  time 
the  space  environment  is  empty,  but  at  times  of  geomagnetic  storms,  plasma  from  the 
sun  interacts  with  the  earth’s  magnetic  field.  In  their  "corpuscular  theory"  this 
interaction  produces  a  magnetic  cavity  in  the  solar  plasma  that  confines  the  earth’s 
magnetic  field  and  produces  magnetic  storms  by  the  establishment  of  a  temporary 
"ring  current"  around  the  earth  in  the  solar  plasma.  We  now  know  that  solar  plasma 
is  always  present  as  the  solar  wind  and  that  their  proposed  magnetic  cavity 
continuously  exists  as  the  magnetosphere.  We  also  know  that  geomagnetic  storms 
and  substorms  are  caused  by  changes  in  the  solar  wind  and  its  interaction  with  the 
magnetosphere. 

Geomagnetic  storms  have  been  recognized  since  the  turn  of  the  century  as 
events  that  impact  our  world  on  a  global  scale.  They  affect  electrical  power 
transmission  systems,  navigational  systems,  radars,  and  cable,  radio,  and  sateUite 
communications.  In  space  they  can  induce  spacecraft  charging  in  satellite  systems 
and  cause  low-orbiting  satellites  to  be  lost.  Being  able  to  understand,  much  less 
forecast,  such  storms  has  been  an  important  goal  of  geophysical  research  for  decades. 

Over  the  past  25  years,  the  primary  focus  of  geomagnetic  storm  research  has 
been  in  the  area  of  auroral  sut  torms.  The  prevalent  thought  has  been  that  if  one 
can  understand  the  substorm,  this  will  lead  to  a  better  understanding  of  the  storm  as 
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a  whole.  While  some  relationships  have  been  derived  for  relatively  quiet  periods, 
these  correlations  break  down  during  times  of  intense  storms.  This  study  will 
investigate  this  relationship  during  intense  storms  and  try  to  determine,  through  the 
use  of  geomagnetic  indices,  how  the  storm  and  substorm  current  ^tems  are  related, 
how  well  they  correlate,  and  if  the  information  contained  in  the  indices  has  any  new 
application  toward  storm  modeling  and  forecasting. 
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CHAPTER  2 

GEOMAGNETIC  STORM  AND  SUBSTORM  PROCESSES 
Magnetospheric  Structure 

As  far  as  storms  and  substorms  are  concerned,  the  magnetosphere  is  where 
most  of  the  action  takes  place.  Understanding  of  the  magnetosphere  and  its 
structure  is  fundamental  to  the  understanding  of  storms  and  substorms. 

According  to  McPherron  [1991],  the  magnetosphere  is  a  composite  of  closed 
and  open  magnetospheric  models.  In  the  closed  model,  viscous  interaction  with  the 
solar  wind  produces  the  magnetospheric  convection  pattern  shown  in  Figure  1  and 
the  corresponding  ionospheric  convection  in  Figure  2  [4xfoni,  1964],  In  the  open 
model,  first  proposed  hyDungey  [1961],  the  IMF  (interplanetary  magnetic  field)  and 
the  earth’s  magnetic  field  merge  at  x-type  neutral  lines  formed  when  the  IMF  is 
southward  (antipmallel).  The  field  lines  connect  at  the  front,  are  carried 
downstream,  and  reconnect  in  the  distant  magnetotail.  This  produces  the  same 
magnetospheric  and  ionospheric  convection  patterns. 

Viscous  interaction  is  always  present  to  some  degree  since  magnetospheric 
convection  is  observed  regardless  of  IMF  orientation.  Open  lines  are  apparently 
always  present  over  the  poles,  since  the  polar  cap  is  always  open  to  the  solar  wind 
[McDiarmid  et  al,  1980].  This  leads  to  the  topology  of  open  and  closed  field  lines 
shown  in  Figure  3.  Open  field  lines  are  connected  at  the  subsolar  point  and  carried 
over  the  poles.  Only  in  the  equatorial  plane  do  closed  field  lines  make  contact  with 
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Fig.  1.  The  magnetospheric  convection  pattern  induced  in  the  equatorial  plane  by 
viscous  interaction  [Axford,  1964]. 


1200 


Fig.  2.  Ionospheric  convection  produced  viscous  interaction  [Axford,  1964]. 
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the  solar  wind,  where  viscous  interaction  drags  the  ISeld  lines  tailward  along  the 
flanks. 

In  Figure  4,  the  composite  motion  of  field  lines  is  depicted  [in  (b)  and  (c) 
motion  due  to  viscous  interaction  is  shaded].  Qosed  field  lines  are  dragged  along  the 
sides  (a)  and  return  just  inside  the  tailward  moving  lines  (b).  The  taihvard-flowing 
field  lines  on  the  flanks  of  the  magnetosphere  form  a  boundary  layer  called  the  low- 


Fig.  3.  Topological  regions  of  the  composite  magnetosphere.  The  solar  wind 
transports  open  field  lines  over  the  poles  (top)  while  viscous  interaction  moves  field 
lines  along  the  flanks  (bottom)  [McPhenvn,  1991]. 
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latitude  boundary  layer  (LLBL).  Open  field  lines  that  merge  at  the  subsolar  point 
are  carried  over  the  poles  (a)  and  flow  earthward  in  the  center  of  the  plasma  sheet 

(b) ,  passing  around  the  plasmapause  and  inside  the  returning  viscous  field  lines.  Part 

(c)  shows  the  corresponding  composite  ionospheric  convection.  Note  that  this  is  done 
for  southward  IMF  -  for  northward  IMF  some  debate  remains  over  whether 
recoimection  occurs  and  over  the  ionospheric  convection  pattern,  although  some 
evidence  shows  that  sunward  convection  occurs  over  the  poles  [Stem  and  Alexeev^ 
1988].  Figure  5  shows  a  cross  section  of  the  magnetotail,  with  the  various  boundary 
layers  depicted.  The  plasma  sheet  actually  has  two  distinct  layers  -  the  Central 
Plasma  Sheet  (CPS)  and  the  Plasma  Sheet  Boundary  Layer  (PSBL).  Because  of 
antiparallel  lines  in  the  two  tail  lobes,  a  current  sheet  (the  cross  tail  current)  flows 
from  dawn  to  dusk  across  the  plasma  sheet  [Ness,  1965]. 

Through  a  combination  of  E  field  drift  and  collisions.  Hall  currents  are 
produced  in  the  ionosphere.  Both  electrons  and  ions  in  the  ionosphere  drift 
according  to  v=ExBlB\  Since  neutrals  are  also  present  and  the  ions  experience 
more  collisions,  this  creates  a  flow  of  electrons  in  the  ExB  direction  and  thus  a 
current  in  the  -ExB  direction.  The  Hall  current  has  the  same  pattern  as  the 
ionospheric  flow  shown  in  Figure  2,  but  in  the  opposite  direction.  The  eastward  and 
westward  auroral  electrojets  are  simply  concentrations  of  these  currents  in  the  high 
conductivity  channels  on  the  night  side  produced  by  auroral  particle  precipitation 
[Akasofu,  1977]. 
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Fig.  4.  Composite  motion  of  magnetic  field  lines  resulting  from  viscous  interaction 
(shaded  region)  and  magnetic  reconnection  {Cowley ^  1982]. 


Charged  particles  trapped  in  the  plasmasphere  (the  area  within  the 
plasmapause)  experience  drifts  due  to  curvature  of,  and  gradients  in,  the  earth’s 
magnetic  field.  The  drifts  are  charge  dependent,  with  positive  partides  drifting 
westward  and  negative  particles  drifting  eastward.  This  produces  a  net  westward 
current  (Figure  6)  that  is  called  the  ring  current  since  regions  of  constant  current 
form  toroids  around  the  earth  [McPherron,  1991]. 
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Fig.  5.  Cross  section  of  the  magnetotail,  showing  the  various  boundary  layers 
[McPherron,  1991]. 


Drift  +  Ion 

+  Ion 


Fig.  6.  A  perspective  view  of  the  dayside  piasmasphere  showing  creation  of  a 
westward  ring  current  by  particle  drifts.  Contours  of  constant  current  form  nested 
toroids  around  the  earth  [McPherron^  1991]. 
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Storm/Substorm  Morphology  and  Relationships 


Definitions 

The  definitions  of  storms  and  substorms  have  evolved  greatly  over  the  years. 
Depending  on  which  author  is  referred  to,  definitions  and  especially  relationships 
between  storms  and  substorms  can  vary  greatly.  For  the  purposes  of  my  research, 
I  will  use  the  definitions  given  by  McPherron  [1991]: 

1)  A  polar  magnetic  substorm  is  a  description  of  the  temporal  and  spatial 
development  of  magnetic  disturbances  during  a  substorm  (this  is  normally  described 
in  terms  of  the  ionospheric  onrent  ^tems).  An  auroral  substorm  is  a  description 
of  the  sequence  of  changes  in  the  aurora  that  accompany  the  polar  magnetic 
substorm.  A  magnetospheric  substorm  is  the  process  that  links  the  magnetosphere 
to  the  ionosphere  and  produces  the  observed  auroral  and  magnetic  activity.  The 
bottom  line  is  that  a  substorm  in  general  is  the  enhancement  of  auroral  currents  and 
particle  precipitation  due  to  a  disturbance  in,  and  increased  energy  input  from,  the 
magnetosphere.  This  disturbance  is  in  turn  rooted  in  changes  in  the  solar  wind’s 
input  into  the  magnetosphere.  Substorms  generally  last  one  to  two  hours. 

2)  A  geomagnetic  storm  is  the  manifestation  in  the  earth’s  surface  magnetic 
field  of  first  a  compression  of  the  magnetosphere  by  enhanced  solar  wind  and  then 
an  intensification  of  the  ring  current  A  geomagnetic  (or  simply  magnetic)  storm  can 
last  for  several  days  and  involves  frequent  substorm  activity. 


10 


Substorm  Phases 

As  a  chiefly  auroral  phenomenon,  substorms  are  best  observed  at  high 
latitudes.  There  is  general  agreement  that  substorms  have  three  identifiable  phases 
as  outlined  below: 

1)  The  growth  phase  is  the  period  during  which  disturbances  in  a  variety  of 
phenomena  become  increasingly  apparent.  There  is  a  gradual  increase  in  the 
eastward  and  westward  electrojets,  and  local  auroral  brightenings  may  occur 
[McPherron,  1979].  In  the  magnetosphere,  the  magnetic  field  strength  in  the  tail 
lobes  increases  and  the  plasma  sheet  thins  in  the  north-south  direction  [Tascione, 
1988]. 

2)  Figure  7a  shows  the  auroral  regions  prior  to  substorm  onset  (quiet 
conditions).  Discrete  auroral  arcs  are  present  in  the  nighttime  auroral  oval  and  in 
the  polar  cap.  At  the  onset  of  the  expansive  phase  the  most  equatorward  auroral  arc 
in  the  midnight  sector  brightens  (Figure  7b).  This  arc  expands  poleward  (the 
poleward  bulge.  Figure  7c)  with  the  region  behind  filled  with  wavelike  auroral  forms. 
The  duskward  edge  of  this  bulge  forms  a  wavelike  disturbance  called  the  westward 
traveling  surge  (Figure  7d).  On  the  dawn  side  the  aurora  forms  wavelike,  omega¬ 
shaped  bands  (omega  bands)  that  drift  eastward  and  the  auroral  arcs  break  up  into 
eastward  drifting,  pulsating  patches  {Akasofii^  1977].  At  the  same  time  as  the  first 
brightening,  there  is  a  strong  increase  in  the  westward  electrojet  The  enhanced 
electrojet  region  expands  westward  with  the  westward  traveling  surge.  The  eastward 
electrojet  also  increases  in  intensity.  At  the  beginning  of  the  expansion  phase  in  the 
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Fig.  7.  Schematic  showing  the  growth  and  decay  of  an  aiux)ral  substorm  \AkasojUy 
1964]. 


magnetotail,  the  tail  field  decreases  and  there  is  an  explosive  acceleration  of  plasma 
towards  the  earth  as  a  portion  of  the  tail  current  is  diverted  along  field  lines  (by  some 
unknown  process).  This  field-aligned  current  connects  through  the  ionosphere  as  the 
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enhanced  westward  electrojet;  collectively  this  system  is  called  the  substorm  current 
wedge  (Figure  8)  [McPherron,  1979]. 

3)  The  recovery  phase  begins  when  poleward  expansion  stops  and  auroral  arcs 
reform  and  drift  equatorward  (Figures  7e  and  f).  The  aurora,  electrojets,  and  the 
magnetotail  field  return  to  their  prestorm  configurations  [McPherron,  1979]. 

Many  substorms  exhibit  much  more  complex  behavior  than  is  described  in  the 
ideal  case.  Multiple  onsets  or  intensifications  in  the  expansion  phase  can  occur  in 
which  many  sudden  intensifications  in  the  electrojet  are  measured.  Pseudo  break-ups 
can  occur,  where  brightenings  of  irregular  auroral  waves  occur  and  expand  poleward. 
These  pseudo  break-ups  last  only  a  few  minutes  and  are  very  weak  and  localized 
[Tascione,  1988]. 


Fig.  8.  A  perspective  view  of  the  substorm  current  wedge  that  develops  during  the 
substorm  expansion  phase  [Ciauer  and  McPherron,  1974]. 
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Storm  Phases 

Magnetic  storms  are  best  observed  at  mid  and  low  latitudes  where  the  primary 
effect  is  a  reduction  of  the  horizontal  magnetic  field  Magnetic  storms  usually  have 
well  defined  phases  as  follows: 

1)  The  initial  phase  begins  with  the  sudden  storm  commencement  (SSC).  The 
SSC  appears  as  a  sudden  worldwide  increase  in  the  horizontal  magnetic  field,  which 
is  due  to  a  compression  of  the  magnetosphere  by  a  shockwave  in  the  solar  wind. 
This  phase  lasts  from  1  to  16  hours  during  which  time  the  horizontal  field  remains 
increased  [Knecht  and  Shuman,  1985]. 

2)  The  main  phase  begins  as  the  horizontal  field  rapidly  decreases  to  below 
pre-storm  levels.  This  is  due  to  the  enhancement  of  the  westward  ring  current  by 
particles  injected  from  the  solar  wind.  Coincident  with  the  decrease,  the  auroral  oval 
expands  equatorward  and  frequent  substorms  occur.  The  main  phase  can  last  from 
several  hours  to  a  few  days  {Knecht  and  Shuman,  1985]. 

3)  The  recovery  phase  consists  of  the  cessation  of  substorm  activity  and  a  quiet 
return  of  the  horizontal  field  to  prestorm  levels.  The  recovery  phase  can  last  for 
several  days  [McPherron,  1991]. 

Storm/Substonn  Relationships 

As  alluded  to  earlier,  there  is  some  disagreement  about  how  storms  and 
substorms  are  related.  Some  believe,  as  Akasofu  [1977,  p.  294]  states,  "...  a 
geomagnetic  storm  period  can  be  identified  as  the  period  when  intense  substorms 
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occur  frequently."  In  1978  C.  T.  Russell  said  (as  quoted  by  Feldstein  [1992,  p.  85]), 
"A  magnetic  storm  is  not  composed  of  successive  substorms.  A  storm  is  something 
different.  Substorms  are  only  incidental,  not  fujidamental,  of  the  storm  development." 

It  seems,  however,  the  majority  believe  that  the  same  mechanisms  that 
produce  substorms  also  produce  storms,  as  is  implied  by  most  authors  of  substorm 
models.  McPherron  [1991]  reasoned  that  whatever  process  that  produces  substorms 
also  creates  enhanced  plasma  convection  in  the  magnetotail.  When  convection  is 
enhanced,  this  increases  the  magnetotail  convection  electric  field  and  this  in  turn 
compresses  the  plasmapause  earthward.  Convecting  particles  now  penetrate  closer 
to  the  earth,  gain  energy,  and  join  the  ring  current.  As  the  tail  electric  field 
fluctuates,  and  the  plasmapause  expands,  the  new  particles  in  this  region  become 
trapped  in  the  plasmasphere  and  thereby  energize  the  ring  current.  This  view  is 
supported  by  other  authors,  such  as  Lyons  and  Schultz  [1989]  and  Fairfield  [1992]. 

Others  believe  the  ionosphere  is  the  primary  source  of  ring  current 
energization  during  storms.  Cladis  and  Francis  [1985]  used  computer  simulation  to 
show  that  storm-time  ring  current  population  can  mainly  be  accounted  for  by  ions 
streaming  out  of  the  polar  ionosphere  along  magnetic  field  lir^s.  Wrenn  [1989] 
studied  particle  lifetimes  in  the  ring  current  to  show  that  a  large  fraction  of  the  ring 
current  is  of  ionospheric  origin. 
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Substorm  Processes 


Driven  Model 

The  driven  model  was  introduced  hyAkasofu  [1981a]  to  show  that  substorms 
are  created  totally  by  a  driven  process.  A  driven  process  is  one  in  which  energy  input 
from  the  solar  wind  is  continuously  fed  into  the  magnetosphere  and  ionosphere  so 
that  no  energy  is  stored  and  then  released  (a  loading-unloading  system).  The  solar 
wind-magnetosphere  dynamo,  which  produces  the  magnetospheric-ionospheric 
convection  and  current  systems,  is  the  mechanism  that  dissipates  energy  from  the 
solar  wind;  however,  detailed  physical  processes  producing  substorms  are  not 
specified. 

Akasofu  defined  the  epsilon  parameter  €  =VBH\m\eri),  where  V is  the  solar 
wind  speed,  B  is  the  magnitude  of  the  IMF,  is  a  constant  proportional  to  the  area 
of  the  magnetopause  through  which  energy  is  input,  and  0  is  the  angle  the  IMF 
makes  with  the  horizontal.  The  «  parameter  represents  the  amount  of  energy  flux 
from  the  solar  wind  into  the  magnetosphere  and  was  derived  from  the  fact  that 
geomagnetic  activity  is  controlled  by  the  strength  and  orientation  of  the  IMF. 
Akasofu  showed  that  c  correlated  well  with  the  total  energy  dissipated  in  the 
magnetosphere  during  substorms,  which  implies  that  substorms  are  a  totalty  driven 
process.  According  to  this  model,  substorms  will  occur  whenever  e  reaches  a  critical 
value  (-10^*  ergs/sec)  and  they  are  a  result  of  increased  power  to  the  solar  wind- 
magnetosphere  dynamo. 
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Near-Earth  Neutral  Line  Model 

The  near-earth  neutral  line  (NENL)  model  is  the  most  widely  accepted  and 
probably  the  most  extensively  developed  model  of  substorm  processes.  The  basic 
premise  of  the  model,  &st  developed  in  the  1970’s  [e.g.,  Russell  and  McPherron,  1972; 
Schindler,  1974;  Hones,  1977]  and  most  recently  reviewed  hy  McPherron  [1991],  is  the 
formation  of  a  second  neutral  line  in  the  magnetotail. 

The  NENL  model  growth  phase  is  due  to  a  driven  process  and  begins  with 
erosion  of  the  magnetopause.  When  magnetic  flux  from  the  solar  wind  increases  and 
exceeds  the  flux  transported  to  the  subsolar  neutral  line  in  the  magnetosphere,  the 
neutral  line  must  move  earthward  to  balance.  The  result  of  this  is  "erosion"  of  the 
dayside  magnetopause  as  successive  layers  of  previously  closed  field  lines  are  opened 
to  the  solar  wind.  These  field  lines  are  transported  over  the  poles  and  flux 
accumulates  in  the  tail  lobes  since  tail  reconnection  does  not  occur  immediately  (this 
delay  is  mostly  due  to  the  time  it  takes  the  solar  wind  to  transport  the  new 
reconnection  discontinuity  to  the  distant  neutral  line).  The  added  flux  causes  the 
distant  tail  lobes  to  grow  in  size  and  causes  an  increase  in  the  intensity  of  the  near- 
earth  lobe  field.  In  response,  the  near-earth  tafl  current  flowing  between  the  lobes 
must  increase  as  the  lobe  field  increases.  Also,  the  near-earth  current  sheet  thins 
because  of  increased  magnetic  pressure  in  the  lobes  and  the  loss  of  magnetic  flux  to 
dayside  reconnection.  Many  of  the  aspects  of  the  growth  phase  are  shown 
schematically  in  Figure  9. 
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The  convection  electric  field  increases  due  to:  1)  a  rarefaction  wave  generated 
by  dayside  reconnection  that  propagates  around  the  earth,  where  in  the  plasma  sheet 
it  accelerates  plasma  sunward;  2)  increased  reconnection  that  allows  the  soleff  wind 
electric  field  to  be  mapped  to  the  polar  cap;  and  3)  the  increase  in  the  near-earth  tail 
magnetic  field.  This  increased  electric  field  moves  the  plasmapause  earthward, 
allowing  the  tail  current  to  move  earthward  and  letting  more  particles  enter  the  ring 
current.  The  ionospheric  electrojets  increase  due  to  enhanced  conductivity  (from  the 
earthward  tail  current  movement  which  increases  particle  precipitation)  and  the 


Fig.  9.  Changes  in  the  magnetosphere  caused  by  a  southward  turning  of  the  IMF 
[McPherron^  1991]. 
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enhanced  electric  field.  This  in  turn  causes  the  Region  1  and  2  currents  to  increase 
by  current  continuity. 

The  basis  of  the  NENL  model  is  the  formation  of  a  neutral  line  in  the 
magnetotail  region  of  10-30  (Figure  10).  Fast  compression  and  rarefaction  waves, 
which  communicate  magnetospheric  changes  to  the  plasma  sheet,  can  reduce  the 
normal  magnetic  field  component  to  zero  in  the  closed  field  line  region  of  the  near- 
earth  plasma  sheet  and  form  an  x-type  neutral  line  [Cornniti,  1985].  Also,  recent 
computer  simulations  using  MHD  processes  support  the  formation  of  a  near-earth 
neutral  line  [Voigt  and  Wolf,  1988;  Bim  and  Hesse,  1991;  Scholer  and  HaiUz,  1991]. 
The  neutral  line  first  forms  in  the  CPS  and  reconnection  begins  there.  Reconnection 
is  initially  slow  since  the  reconnection  rate  is  proportional  to  the  Alv6n  velocity,  which 
is  slow  in  the  CPS.  As  reconnection  continues,  it  reaches  the  PSBL  and  the  rate 
increases  due  to  the  increased  Alv6n  velocity.  Eventually  the  last  closed  field  line 
that  defines  the  plasma  sheet  boundary  is  severed  and  the  recoimection  rate 
dramatically  increases  in  the  tail  lobes.  In  the  tail  anti-earthward  of  the  near-earth 
neutral  line,  a  plasmoid  (a  'bubble"  of  closed  magnetic  field  lines)  is  formed  (Figure 
11).  The  onset  of  the  expansion  phase,  which  is  due  to  a  loading-unloading  process 
(the  accumulation  and  release  of  lobe  magnetic  flux),  is  defined  as  the  moment  the 
last  closed  field  line  is  cut  and  the  plasmoid  begins  to  retreat.  The  plasmoid  retreats 
anti-sunward,  as  shown  in  Figure  11,  because  of  the  tension  in  the  IMF  field  lines 
that  now  separate  it  fi’om  the  near-earth  neutral  line.  The  NENL  model  is  largely 
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based  on  distant  tail  observations  that  seem  to  show  the  passage  of  a  plasmoid 
[Hones  et  al,  1984;  Baker  et  al,  1987;  Nishida  el  al,  1988;  Slavin  et  aL,  1989]. 

The  formation  of  a  plasmoid  also  requires  some  kind  of  current  diversion  from 
the  reconnection  region.  To  produce  an  x-line  geometry,  the  tail  current  in  the  x-line 
region  must  be  reduced.  This  can  be  done  by  introducing  field-aligned  currents  that 
divert  current  away  from  the  x-line  region  and  close  through  the  ionosphere  (Figure 
12).  This  is  how  the  substorm  current  wedge  is  formed. 

Eventually  the  NENL  retreats  tailward,  as  depicted  at  the  bottom  of  Figure 
11,  during  the  recovery  phase.  There  is  no  clear  cause  of  this  retreat,  but  when  the 
IMF  turns  northward  and  reconnection  ceases,  there  is  flux  in  the  magnetotail  waiting 
to  be  reconnected.  The  simplest  way  to  solve  this  problem  is  for  the  neutral  line  to 
move  tailward,  adding  excess  flux  to  the  plasma  sheet  behind  it. 


Fig.  10.  The  formation  of  a  near-earth  neutral  line  within  the  CPS  during  the  late 
growth  phase  [McPherron,  1991]. 
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Electrojet  Field  Near 


Fig.  12.  The  diversion  of  the  tail  current  through  the  ionosphere  during  the  substorm 
expansion  phase  [McPhenon^  1991]. 
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CHAPTERS 
DATA  DESCRIPTION 


Geomagnetic  Indices 


Dst  Index 

The  Dst  index  was  first  derived  both  by  Sugiura  [1964]  and  Kertz  [1964]  and 
was  conceived  as  a  measure  of  the  intensi^  of  the  equatorial  ring  current.  Since  the 
ring  current  is  assumed  to  be  mostfy  a  symmetric,  equatorial,  zonal  current,  its 
magnetic  field  at  the  surface  is  parallel  to  the  dipole  axis  and  seen  mainly  in  the 
horizontal  component  (JT)  of  the  magnetic  field  [Rangarajan,  1989].  The  method  of 
Sugiura  is  essentially  that  which  is  in  use  today,  in  which  the  average  disturbance  of 
the  horizontal  field,  D,  is  made  up  of  two  parts:  1)  a  universal-time  part  Dst;  and  2) 
a  local-time  dependent  part  DS.  For  a  uniform  longitudinal  distribution  of  stations, 
the  mean  value  of  D  at  any  instant  will  not  be  contaminated  by  DS,  so  that  essentially 
Dst  is  measured.  And  since  the  spatial  and  temporal  variations  in  ring  current 
intensity  are  not  rapid,  a  dense  network  of  stations  is  not  needed.  Currently,  Dst  is 
measured  by  four  stations  (Table  1)  that  are  located  far  enough  from  the  magnetic 
equator  so  that  they  are  not  affected  by  the  equatorial  electrojet,  yet  close  enough 
to  the  equator  to  minimize  auroral  effects  [Rangarajan,  1989]. 

In  computing  the  index,  hourly  observations  at  each  station  are  taken.  The 
diurnal  solar-quiet  (Sq)  variation  is  computed  at  each  station  for  each  month  for  the 
5  quietest  days,  then  from  the  series  of  monthly  Sq  for  the  year,  a  double  Fourier 
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Table  1.  Coordinates  of  the  magnetic  stations  of  the  Dst  network  [adapted  from 
Rangarajan,  1989]. 


Station 

Geographic 

Latitude 

Geographic 

Longitnde 

Dipole 

Latitude 

Dipole 

Longitude 

Honolulu 

21.32 

-158.00 

21.36 

267.97 

San  Juan 

18.12 

-66.16 

29.41 

4.65 

Hermanns 

-34.42 

19.23 

-33.59 

81.99 

Kakioka 

36.23 

140.18 

26.31 

207.25 

series  is  expanded: 

EE'*.'  cos(wr+aJ  cos (nM+/J^  (1) 

*•1 


where  T  is  local  time  and  M  is  the  month.  One  then  has  12x24  experimental  data 
to  determine  the  48  unknown  coefficients.  The  hour  values  of  a  synthetic  Sq 
variation  are  computed  from  the  Fourier  series  and  are  subtracted  from  the  original 
hourly  values.  The  reference  level  at  each  station  is  determined  a  parabola  fitted 
to  the  annual  means  of  H  for  the  5  quietest  days  each  month  {Mayaud,  1980]. 

Although  a  synthetic  Sq  is  used  to  remove  regular  variations,  some  residual 
variation  may  remain  in  the  index.  Also,  the  effects  of  the  partial  ring  current  and 
magnetospheric  compression  by  the  solar  wind  are  not  removed.  This  leaves  the 
symmetric  ring  current,  partial  ring  current,  and  magnetospheric  compression  as  the 
primary  sources  for  Dst.  During  quiet  periods  magnetospheric  compression  is 
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removed  along  with  Sq,  but  during  storms  the  effect  can  be  10  to  40  nT 
[Bawnjohann,  1986]. 

Auroral  Indices 

The  auroral  indices  AE,  AL,  AU,  and  AO  (hereafter  called  the  AE  indices) 
were  first  introduced  by  Davis  and  Supura  [1966]  to  meastire  the  auroral  electrojets 
and  are  discussed  in  detail  by  Mayaud  [1980].  Using  a  chain  of  12  auroral  stations 
(Figure  13),  a  quiet  time  reference  level  is  derived  for  each  month  at  each  station 
from  the  5  quietest  days.  This  level  is  subtracted  from  the  2.5  or  1-minute  H  data 
for  each  station,  the  data  is  superposed  in  universal  time,  and  the  extreme  positive 
and  negative  values  at  each  instant  define  AU  and  AL,  respectively.  The  2.5  or  1- 
minute  values  are  averaged  to  produce  hourly  values  which  are  the  ones  used  in  this 
study.  AL  should  measure  the  intensity  of  the  westward  electrojet,  while  yll/  should 
measure  the  intensity  of  the  eastward  electrojet.  However,  these  indices  also  contain 
contributions  from  any  other  zonal  currents  (in  the  ionosphere  and  beyond),  mainly 
the  ring  current,  so  AE  is  defined  as  AU—AL  so  as  to  remove  any  symmetric  zonal 
contribution.  AO,  defined  as  (AU+AL)I2,  is  then  intended  to  be  an  approximate 
measure  of  the  equivalent  zonal  current.  Since  there  is  asymmetry  between  the 
eastward  and  westward  electrojets,  however,  .<40  really  measures  this  asymmetry  plus 
zonal  current  effects. 

Some  problems  with  the  AE  indices  that  must  be  remembered  are  [Mayaud, 


1980]: 
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Fig.  13.  The  chain  of  AE  stations.  Geographic  coordinates  are  indicated  by  solid 
lines,  and  geomagnetic  coordinates  are  shown  by  plus  signs  [from  Data  Book  21^ 
WDC-C2  for  Geomagnetism,  Kyoto,  Japan,  1992]. 


1)  Underground  effects  induced  by  the  rapidly  fluctuating  currents. 

2)  The  direction  of  H  (magnetic  north)  in  which  the  measurements  are  made 
may  not  be  normal  to  the  electrojets. 

3)  The  present  network  of  stations  is  not  ideal.  Since  longitudinal  gaps  are 
present,  a  low  AB  does  not  preclude  substorm  activity. 
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4)  Ring  current  effects  are  present  mAV  andylL  and  are  important  when  Dst 
is  large. 

5)  The  daily  regular  variation  of  H  is  not  eliminated.  Although  this  is  not 
significant  during  active  periods  (the  daily  variation  is  tens  of  nT  while  the  electrojet 
variations  are  hundreds  of  nT),  the  daily  variation  itself  can  define  .4(7  and  .<4L  during 
quiet  periods. 

6)  The  indices  caimot  distinguish  between  movement  of  the  auroral  oval  away 
fi’om  (or  towards)  the  station  and  weakening  (or  strengthening)  of  the  electrojet. 

7)  According  to  BertheUer  [1976],  the  component  of  the  IMF  makes  some 
contribution  to  the  AE  index. 

Akasofri  etal  [1983]  addressed  the  question  of  the  accuracy  of  Xh.tAE  indices 
by  comparing  them  to  AE  from  a  chain  of  70  stations  used  during  the  IMS 
(International  Magnetospheric  Study).  They  found  that  the  relative  accuracy  of  the 
AE  indices  became  progressively  worse  for  lower  values  (although  the  overall 
correlation  was  -0.93)  and  suggested  a  critical  value  of  250  nT,  below  which  one 
should  be  cautious  of  the  value  of  AE.  Overall,  however,  the  advantages  of  the  AE 
indices  far  outweigh  the  deficiencies  as  the  indices  are  directly  related  to  physical 


processes. 
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Data  Period  and  Data  Sets 


Storm  Criteria  and  Data 

The  problem  of  how  to  define  a  storm  is  not  simple  since  one  can  be  defined 
in  so  many  ways.  They  can  be  classified  in  terms  of  the  indices  Ap,  ap,  Kp,  or  Dst. 
A  storm  can  also  be  defined  in  terms  of  its  efiects,  both  terrestrial  and  orbital.  For 
this  stucty  I  will  adopt  but  modify  the  criteria  used  by  NOAA’s  Space  Environmental 
Services  Center  (SESC).  Their  criterion  for  a  severe  geomagnetic  storm  S&Ap'S:.  100, 
however,  I  will  v&cAp*  ^  100  (Ap  is  taken  fi^om  0000  to  2400  GMT,  while  is  the 
maximum  possible  value  over  any  24  hour  period).  I  will  also  supplement  this  with 
the  additional  criterion  oiDst  s  -200  nT  to  insure  no  storms  with  large  ring  current 
intensifications  are  missed. 

The  data  used  for  most  of  this  study  is  geomagnetic  indices  and  IMF  data 
from  solar  cycle  21,  from  1978  to  1986.  Table  2  shows  data  for  the  21  storms  that 
meet  the  above  criteria,  and  Figure  14  shows  the  superposed  Dst  traces  for  the  21 
storms. 

Storm  Examples 

Figure  15  shows  an  example  oiAE  and  Dst  for  a  "normal"  storm.  In  the  top 
figure  is  the  1-minute  AE  for  the  storm.  It  exhibits  very  rapid  fluctuations  and  has 
an  almost  noisy  appearance.  The  middle  figure  is  the  1-hour  AE  index,  which  is  an 
average  of  the  1-minute  data.  The  trace  still  shows  that  the  AE  index  is  subject  to 
rapid  temporal  changes  since  the  substorm  time  scale  is  relatively  short  (1-2  hours). 


Table  2.  Summary  of  data  for  the  21  storms  selected  from  solar  cycle  21 


Date 

SI  Day 
Fia7 

May  1-4, 1978 

142 

Aug  27-30, 1978 

135 

Sep  27-30, 1978 

152 

Apr  34, 1979 

178 

Apr  24-28, 1979 

177 

Dec  19-22, 1980 

193 

Mar  5-8, 1981 

209 

Apr  11-14, 1981 

212 

Jul  24-27, 1981 

206 

Mar  1-4, 1982 

202 

Jul  13-16, 1982 

172 

Aug  1982 

166 

Sep  5-8,1982 

163 

Sep  21-24, 1982 

163 

Jan  9-12, 1983 

155 

Feb  44, 1983 

128 

Apr  25-29, 1984 

128 

Sep  22-25, 1984 

77 

Nov  15-18, 1984 

73 

Apr  19-23, 1985 

78 

Fa>  7-10, 1986 

76 

Siorm 

Nninbcr 


Max 

ShrKp 


Fig.  14.  Superposed  Dst  traces  for  the  21  storms  from  solar  cycle  21.  The  peak  of 
each  trace  occurs  at  t=0. 
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At  the  bottom  is  a  classic  Dst  storm  profile,  where  each  storm  phase  is  visible.  From 
an  initial  baseline  of  low  activity  (Osr-O),  the  trace  exhibits  a  sharp  increase 
associated  with  magnetospheric  compression  (the  sudden  storm  commencement). 
The  field  then  decreases  rapidly  as  the  ring  current  is  intensified  (main  phase), 
reaches  its  maximum,  and  decays  slowly  back  to  the  baseline  level  (recovery  phase). 
It  is  evident  firom  Figure  15  that  very  different  time  scales  are  at  work  in  the  two 


current  systems. 
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CHAPTER  4 

COMPARISONS  OF  DST  AND  AE 

Rationale 

While  a  relationship  between  the  auroral  electrojets  and  the  ring  current 
seems  plausible,  since  both  are  driven  by  the  same  basic  source  (the  solar  wind),  very 
rarely  has  this  relationship  been  explored.  The  veiy  nature  of  what  the  relationship 
might  be  is  unknown,  and  some  doubt  one  even  exists.  Are  substorms  a  fundamental 
part  of  storm  development,  or  are  they  coincident  but  not  necessary?  An 
investigation  of  the  association  between  the  storm  and  substorm  current  systems 
should  be  able  to  provide  at  least  a  start  on  answering  this  question. 

Previous  Work 

One  of  the  first  direct  comparisons  of  Dst  and  AE  was  made  by  Davis  and 
Parthasarathy  [1967].  They  found  that  1)  the  onset  of  AE  activity  precedes  DR 
(similar  to  Dst)  enhancement  by  as  much  as  15  hours  for  geomagnetic  storms;  2)  the 
amplitude  of  the  maximum  DR  is  directly  proportional  to  the  sum  of  MatAE  values 
for  the  previous  10  hours  (th^  achieved  a  correlation  coefficient  of  0.82  for  32 
magnetic  storms  in  1958);  and  3)  the  energy  injection  function  to  the  ring  current 
(derived  by  them)  looks  very  similar  to  the  time  variation  of  AE,  indicating  that  both 
are  energized  by  the  same  process. 
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Akasofu  [1981c]  made  a  direct  comparison  of  the  AE  and  Dst  indices  for 
several  different  solar  rotation  periods.  Figure  16  shows  a  scatter  plot  of  a  quiet 
rotation  period  and  Figure  17  shows  a  plot  for  a  period  that  contained  a  geomagnetic 
storm.  Akasofu  concluded  that  while  there  is  a  linear  relationship  at  relatively  quiet 
periods  (Dst  >  -50  nT),  during  storms  there  is  no  longer  a  relationship  between  the 
two  indices  since  during  the  lowest  Dst  values  the  AE  index  could  be  relatively  low. 

Siscoe  [1982]  used  an  electrical  circuit  analogy  to  study  the  energy  coupling 
between  the  ring  current,  ionospheric  currents,  and  field-aligned  currents.  He  found 
that  while  AE  varies  linearly  with  the  component  of  the  IMF,  Dst  varies  as  5/,  so 
that  the  two  indices  would  necessarily  have  different  responses  to  a  southward  IMF. 


FEB.  10  .  1974  -  MAR. 8  ,1974 


DST 

COMMAS 


Fig.  16.  The  relationship  between  and  Dst  indices  for  the  Bartel’s  rotation  period 
1922  (February  10-March  8,  1974)  [Akasofu,  1981b]. 
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Fig.  17.  The  relationship  between..4£  andZlsr  indices  for  the  Bartel’s  rotation  period 
1910  (March  23-April  18, 1973)  \Akasofu,  1981b]. 


While  not  using  the  AE  index,  Wrenn  [1989]  compared  3-hour  ap  to  Dst  by 
using  a  weighting  factor  r  applied  to  ap.  This  attenuation  factor  r  was  first 
introduced  by  Wrenn  [1987]  in  comparing  ap  to  ionospheric  parameters.  This  factor 
(which  can  vary  between  0  and  1)  added  increased  weighting  to  previous  values  and 
imposed  an  exponential  decay  time  constant  of  3/(1-t)  hours.  The  new  time-weighted 
index  was  given  by: 

Although  Wrenn  compared  ap(T)  to  DsU  he  disregarded  intense  storms  eliminating 
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data  100  hours  after  Dst  <  ISO  nT  and  6  hours  after  a  sudden  storm  commencement 
From  his  comparisons  he  concluded  that  ring  current  particles  have  longer  lifetimes 
during  solar  maximum  than  solar  minimum  and  that  a  significant  portion  of  ring 
current  particles  is  of  ionospheric  origin. 

In  the  recent  review  of  current  storm  and  substorm  research  by  Fairfield 
[1992],  he  gave  Figure  18  as  a  summary  of  the  current  thinking  on  substorm  research. 
This  figure  indicates  that  AE  is  dominated  by  directly  driven  processes  (i.e., 
magnetospheric  convection)  while  the  actual  impulsive  substorm  component  (from 
loading-unloading  processes)  is  superposed  on  the  directly  driven  part  This  view  is 
shared  by  Lee  et  aL  [1985]  and  Uu  et  at  [1988]  who  both  state  that  the  AE  index  is 
composed  of  an  enhanced  level  component  that  is  directly  driven,  with  an  impulsive 
component  (fi-om  plasmoid  formation)  superposed  on  this  background  level  Not 
shown  in  Figure  18  is  the  ring  current,  which  is  largely  believed  to  also  be  directly 
driven  by  magnetospheric  convection.  This  implies  that  some  relationship  should 
exist  between  Dst  and  the  directly  driven  component  olAE. 

Smoothing/Weighting  Techniques  Applied  to  AE 

From  Akasofu’s  work,  it  seems  clear  that  no  elementary  relationship  can  exist 
between  AE  and  Dst  because  the  ciurent  systems  they  measure  have  fundamentally 
different  time  constants.  While  Dst  exhibits  a  fairly  fast  rise  and  very  slow  decline, 
AE  can  fluctuate  greatly  over  a  period  of  hours  (see  Figure  15),  hence  the  large 
scatter  in  Figure  17.  If  one  were  able  to  compensate  for  the  different  time  constants. 
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Fig.  18.  Summary  of  much  of  current  thinldng  on  magnetospheric  storms  [Fairfield, 
1992]. 


a  clearer  picture  might  emerge.  This  was  attempted  on  a  limited  basis  by  Davis  and 
Parthasarathy  when  they  found  that  the  maximum  DR  of  a  storm  was  related  to  the 
integrated  AE  activity  over  the  past  10  hours.  In  addition,  by  looking  at  Figure  18, 
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one  might  also  find  a  better  relationship  ilDst  is  compared  to  only  the  directly  driven 
component  of  AE. 

A  simple  space-centered  smoothing  technique  was  applied  to  the  lAxowAE 
data  for  each  of  the  21  storms  selected  for  this  study.  This  procedure  was  used  to 
attempt  to  estimate  the  trend  of  the  data  and  recover  an  average  baseline  that  might 
better  approximate  the  directly  driven  currents.  An  example  of  this  technique  for 
different  smoothing  iterations  is  shown  in  Figure  19. 

The  technique  of  time-weighted  accumulations  applied  by  Wrenn  [1987]  to  ap 
appears  to  be  able  to  resolve,  at  least  partially,  the  problem  of  the  differing  time 
constants.  K  the  t  technique  were  applied  to  AE,  this  might  produce  a  better 
correlation  to  Dst,  especially  during  storms.  The  following  formula  was  used  to  apply 
the  T  weighting  to  the  AE  indices: 

AE,(t)  =  AE,.j(r)+(l-T)[AE,-AE,.i(r)]  (3) 

This  was  done  for  several  values  of  t,  examples  of  which  are  shown  in  Figure  20.  It 
is  apparent  in  the  figure  that  as  the  value  of  r  increases,  more  weight  is  given  to 
previous  indices  and  the  curve  becomes  smoother  and  flatter.  This  has  the  physical 
significance  of  injecting  the  auroral  energy  into  a  system  with  a  much  slower  response 
time  or  time  constant.  The  new  time  constant  for  a  He  decay  can  be  defined  as  « 
1/(1-t)  hours,  which  for  a  t  of  0.9  is  10  hours. 


1  Hr  (nT)  1  Hr  (nT)  1  Hr  AE,,„  (nT)  1  Hr  AE,„,  (nT) 


July  13-16,  1982 


Fig.  20.  Examples  of  time-weighted  accumulations  of  AE  (thin  line)  for  r  values  of 
0.6,  0.7,  0.8,  and  0.9  to  achieve  AE(t)  (bold  line). 
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Before  proceeding,  a  note  should  be  made  on  the  use  of  1-hour  AE  versus  1- 
minute  AE.  Applying  Fourier  analysis  to  the  1-minute  AE  gives  several  interesting 
results  relating  to  the  1-hour  index.  When  a  Fourier  transform  of  the  l-minute  index 
is  reconstructed  with  the  higher  frequencies  removed,  the  data  becomes  very  close 
to  the  1-hour  index.  Therefore,  the  1-hour  index  can  arguably  be  said  to  represent 
the  lower  frequency  variations  of  the  electrojets.  Also,  as  the  high-end  frequency  cut¬ 
off  is  lowered,  the  reconstructed  index  approaches  that  derived  by  simple  space 
centered  smoothing  on  the  1-hour  index,  so  that  this  methcxl  approximates  even 
lower  frequency  components  of  the  electrojet  variations.  When  t  weighting  is  applied 
to  the  l-minute  index  (as  well  as  reconstructed  indices  with  higher  frequencies 
removed),  it  is  almost  identical  to  that  derived  from  the  1-hour  AE  index.  All  of 
these  results  are  discussed  in  more  detail  (with  examples)  in  Appendix  B. 

Correlations  with  Dst 

Applying  the  two  above  techniques  to  the  AE  series  of  indices  produced  two 
new  sets  of  adjusted  indices,  referred  to  asXE^  (for  space  centered  smoothing)  and 
AE(t)  (for  time-weighted  accumulations).  These  new  sets  for  the  21  storms  were 
generated  for  48  hours  before  and  72  hours  after  the  minimmn  Dst  value,  then 
compared  to  Dst.  After  trial  and  error  to  determine  the  best  values  to  use,  the 
comparisons  were  made  for  50  iterations  of  averaging,  and  t  values  of  0.75  and  0.9. 
Plots  of  these  for  each  of  the  storms  are  shown  in  Appendix  A.  Correlation 
coefficients  were  computed  for  each  of  the  indices  AE,  AL,  AU,  and  AO  compared 
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to  Dst,  then  each  index  with  the  averaging  and  t  weighing.  As  can  be  seen  from  the 
plots  in  Appendix  A,  in  some  cases  there  is  an  apparent  time  shift  between  the  two 
curves.  Therefore,  the  curves  were  time  shifted  to  achieve  a  maximum  correlation. 
The  resulting  correlation  coefficients  are  shown  in  Figures  21-24,  each  with  a 
histogram  of  the  time  shifts  needed  to  achieve  maximum  correlation  (a  positive  ofrset 
shifts  the  AE  indices  forward).  The  figures  show  that  ih^AO  set  of  indices  has  the 
best  correlations,  followed  by  AL,  with  the  t  weighing  of  0.9  being  the  best  for  all 
indices.  The  time  shift  plots  demonstrate  that  AE  correlates  with  Dst  using  a  time 
shift  centered  on  -6  hours,  while  a  t  of  0.9  reduces  the  time  shift  to  near  0.  Similar 
results  are  found  for  XL.  The  results  forXt/  are  less  consistent,  with  the  maximum 
correlation  for  the  straight  index  often  shifted  by  24  hours  (the  maximum  allowed), 
while  XO  is  very  consistent  in  being  clustered  near  0  for  most  values.  The  similarity 
of  AE  and  XL  is  expected,  since  XL  makes  a  much  greater  contribution  than  AU. 
The  poorer  results  olAU  can  be  explained  by  the  fact  that  the  eastward  electrojet 
is  less  correlated  to  substorm  activity  and  is  not  well  understood.  The  good 
correlation  of  XO  is  justified  since  AO  is  intended  to  measme  the  zonal  currents,  of 
which  the  ring  current  is  a  part.  XL  should  also  be  slightly  better  thanXE  since  X/ 
contains  a  ring  current  contribution  and  AE  does  not. 

The  most  significant  result  is  the  close  correlation  between  Dst  and  AE(t) 
(ranging  from  0.73  to  0.93  except  for  storm  3,  which  was  0.58)  andXL(T)  (ranging 
from  0.78  to  0.95  except  for  storm  3,  which  was  0.67)  for  t=0.9.  By  introducing  a 
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time  constant  of  ~10  hours  to  the  AE  and  AL  data,  one  can  achieve  very  close 
correlations  with  Dst  during  intense  storms. 
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Fig.  21.  Correlation  coefficients  (a)  and  time  shifts  needed  to  achieve  mairinniin 
correlation  (b)  for  AE,AE^  and.4E(T). 
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Fig.  22.  Correlation  coefficients  (a)  and  time  shifts  needed  to  achieve  maximum 
correlation  (b)  ioxAL,AL^  mdALif). 
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Fig.  24.  Correlation  coefficients  (a)  and  time  shifts  needed  to  achieve  maximum 
correlation  (b)  for  AO,  AO^  ^ndAO(T). 
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CHAPTERS 

DST  MODELING  FROM  AE  AND  AL 
Previous  Modeling  Efforts 

Several  methods  for  modeling  Dst  (and  hence  the  ring  current)  have  been 
developed  over  the  past  30  years.  An  early  effort  by  Kamide  and  Fukushima  [1971] 
attempted  to  use  AE  times  an  exponential  function  to  model  the  energy  input  into 
the  symmetric  ring  current.  By  varying  the  decay  time  in  the  equation  for  energy 
change  of  the  ring  current,  they  were  able  to  represent  Dst  fairly  well.  However,  this 
was  done  for  only  one  storm  and  no  further  work  in  this  area  was  ever  pursued. 

Akasofii  [1981b]  used  his  e  parameter  to  model  AE  and  Dst  to  show  that  both 
are  directly  driven  by  the  solar  wind.  He  achieved  good  agreement  for  a  few  storms 
and  indicated  the  possibility  of  forecasting  these  indices  from  IMF  measurements. 

The  most  current  methods  are  reviewed  by  FeJdstein  [1992].  All  of  these 
methods  are  alike  in  that  they  use  IMF  parameters  to  model  Dst.  Feldstein  tested 
several  of  them  and  determined  that  the  method  developed  by  Pisarsky  et  aL  [1989] 
was  one  of  the  best.  The  Pisarsky  et  al.  method  was  applied  to  those  storms  of  the 
21  which  had  close  to  a  complete  set  of  IMF  data  and  the  resiilts  are  given  in 
Appendix  C.  One  problem  with  the  Pisarsky  et  al.  method  is  that  it  requires  a  prior 
knowledge  of  the  peak  Dst  so  that  it  cannot  model  Dst  in  real  time.  Other  methods, 
such  as  that  introduced  by  Burton  et  aL  [1975],  could  model  Dst  in  real  time  if  IMF 
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were  always  available.  However,  another  difficulty  is  that  IMF  data  is  not  always 
available  for  a  storm,  much  less  available  without  delay. 

General  Relationship  Between  AE/AL  and  Dst 

Given  the  high  correlations  found  between  the  AE  indices  and  Dst  in  the 
previous  chapter,  and  the  close  representations  to  the  Dst  traces  shown  in  ^pendix 
A,  it  is  possible  that  AE  could  be  used  to  model  Dst  and  demonstrate  a  current 
system  relationship,  with  AE(t)  and  AL(t)  (with  t=0.9)  being  the  best  candidates 
fAO(T)  will  not  be  considered,  even  though  it  had  the  best  correlation,  since  it  has 
a  smaller  lead  time  (see  Figure  24)  and  would  be  less  useful  for  forecasting  —  it  also 
has  less  physical  significance].  The  year  1982  was  chosen  for  use  in  determining  the 
relationship  between  the  indices  since  it  contained  a  large  number  of  intense  storms. 
AE(t)  and  AL(t)  with  t=0.9  [afterwards  referred  to  as  AE(90)  and  AL(90)]  were 
then  calculated  for  the  entire  year  and  compared  to  Dst.  The  correlation  coefficients 
for  the  year  were  -0.751  for  A£'(90)  and  0.792  for.4L(90),  compared  to  -0.571  for 
AE  and  0.584  for  AL.  Figure  25  shows  a  scatter  plot  for  the  entire  1982  data  set 
It  is  apparent  in  the  plot  that  a  linear  relationship  exists  for  lower  values,  while  an 
exponential  function  might  better  represent  the  higher  values. 

A  least-squares  fit  was  performed  on  the  data  points  in  Figure  25.  Since  the 
higher-valued,  nonlinear  points  were  relatively  few,  they  did  not  affect  the  fitting. 
The  resulting  equations  were: 
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Dst0iE(9Q))  =  [-0.147  Xi4£(90)]  +  20.2  (4) 

Dst(4U90))  =  [0.224  xAL(90)]  +  17.2  (5) 

Values  of  Dst,  AE(90),  and  AL(90)  were  then  taken  for  65  Dst  peaks  during  the 
year.  These  values  were  plotted  on  a  logarithmic  scale  and  fitted  with  a  least-squares 
fit  to  produce  the  following  relations: 

DstiAE(90))  =  -13.2xltf®““**^>l  (^) 

Dst(4L(90y)  =  -13.8xl0t"®“‘^^~^»  0) 


The  peak  values  used  and  the  least-squares  fit  are  shown  in  Figure  26. 

To  use  these  two  relationships  together,  it  must  be  determined  when  to  use 
each.  For  AE,  the  curves  for  the  two  equations  intersect  at  y4£(90)=450  (Dst=-A6) 
and  AE(90)=575  (Dsr=-64),  while  for  XL  they  intersect  at  XL(90)=-220  (Ds/=-32) 
and  XL(90)=-450  (Dst=-83).  For  both  cases,  the  higher  values  were  chosen  since 
from  Figure  25  the  linear  relation  clearly  holds  up  to  those  points.  This  results  in  the 
following  composite  relationships: 


Dst(AEC90))  =  [-0.147  Xi4£(90)i  +  20.2;  i4£(90)s575 
Dst(AEi90))  =  -13.2x10*°°°“’’^**®^;  X£(90)>575 


Dst(4L(^90))  =  [0.224  Xi4I<90)]  +  17.2;  i4I(90)i-450 

Dst(AlJ(,90))  =  -13.8xl0l'^®®”*’‘^>l;  AL(90)<-450 


(8a) 

(8b) 

(9a) 

(9b) 
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When  equations  (8)  and  (9)  are  applied  to  the  data  for  1982,  the  correlation 
coefficient  for  D5r(AE(90))  is  0.775  and  for  D5t(AL(^))  is  0.793.  Monthly  plots  for 
Dst  versus  Dst{AE(^))  and  Ds/(AL(90))  for  1982  are  in  Figures  27  and  28.  In  these 
plots  it  can  be  seen  that  for  the  majority  of  the  time,  7>y/(AE(90))  and  D5/(/iL(90)) 
give  quite  an  accurate  representation  of  Dst. 

Since  (8)  and  (9)  were  developed  from  the  1982  data,  one  would  expect  a 
reasonable  coirelation.  To  test  the  validity  of  the  derived  relations,  (8)  and  (9)  were 
applied  to  the  other  years  of  this  study  (1978-1986)  and  the  results  are  given  in  Table 
3.  For  most  years  the  relationships  work  just  as  well,  and  for  some  years  (1981  and 
1984)  they  are  even  superior,  with  correlations  reaching  0.8  or  better. 

Relationships  Between  AE/AL  and  Dst  for  Storms  1-21 

As  can  be  seen  in  Figures  27  and  28,  the  general  relationships  (8)  and  (9) 
model  Dst  quite  well  for  most  of  the  year.  However,  dining  intense  storm  periods 
(8)  and  (9)  sometimes  work  well  and  other  times  do  not.  To  try  to  determine  the 
nature  of  these  discrepancies,  the  relationships  between  Dst  and/lE(90)  andAL(90) 
were  investigated  in  detail  for  storms  1-21.  The  functional  relationships  derived  are 
of  linear  or  exponential  form  of  these  two  types: 

Dst  =  [axiindex)]  +  b  (1®) 


Dst  = 


(11) 


ms 


Table  3.  Correlation  coefficients  for  Dst  versus  Dst  as  a  function  of 
AE{9G)  and  AL(90)  for  the  years  1978-1986. 
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Year 

Correlation  Coefficients  | 

Dst(AE(90)) 

Dst(AL(90)) 

1978 

0.7645 

0.8012 

1979 

0.7032 

0.7654 

1980 

0.7113 

0.7920 

1981 

0.7722 

0.8180 

1982 

0.7753 

0.7927 

1983 

0.7224 

0.7755 

1984 

0.8081 

0.8242 

1985 

0.7525 

0.7898 

1986 

0.7453 

0.6050 

The  coefficients  for  these  relations  for  each  storm  are  shown  in  Table  4  and  the 
functions  are  plotted  in  comparison  to  relations  (8)  and  (9)  in  Figure  29.  The 
resulting  correlation  coefficients  are  given  in  Table  5.  Possible  reasons  for  these 
functional  differences  are  discussed  in  the  next  chapter. 
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Table  4.  Functional  relationship  between  Dst  and  AE(90)IAH90)  for  storms  1-21. 


1  Sionn 

FuKtiooal  Dependence  of  Dst  on  AE(90) 

Fanctlonal  Dependence  of  Dst  on  AL(90)  | 

!■ 

eq.  (8a);  AE(90)a575 

eq.  (ll)o4Uic— 19.6,  d-.000891;  AE(90)>575 

eq.(9a);  AL(90)a-450  | 

eq.(ll)witbc— 29.7,d— .00099%AL(90)<-450  | 

2 

eq.  (8a);  AE(90)aS7S 

eq.  (11)  with  c»  -3.47,  d-.0022;  AE(90)>575 

eq.  (9a);  AL(90)2-450  | 

eq.(ll)withc—9.01,d— 30214;  AL(90)<-450  | 

3 

eq.  (8a);  AE(90)a575 

eq.  (11)  with  c<— 1.87,  d-.00267;  AE(90)>575 

eq.(9a);  AL(90)&-450  1 

eq.  (11)  with  c-  -349,  d-  -.00306;  AL(90)<  -450  | 

:■ 

eq.  (10)  with  a— .209, 8-25.1;  AE(90)a550 

eq.  (10)  with  a-  -.60,  b-240;  AE(90)>550 

eq.  (10)  with  a-330, 8-193;  AL^)2-393  | 

eq.  (10)  with  a-.709, 8-169;  AL(90)<  -393 

s 

eq.  (8a);  AE(90)S575 

eq. (11) withe— 9.95, d=.00141;  AE(90)>575 

eq.  (9a);  AL(90)S-450 

eq.(ll)withe— 5.19,  d— .00268;  AL(90)<-450 

6 

eq.  (10)  with  a— —.605,  b-90.7 

eq.  (10)  with  a-387, 8-58.7 

7 

eq.  (10)  with  a— 336,  b-16a 

eq.  (10)  with  a«.434, 8-17.4 

8 

eq.  (10)  with  a— 311,  b-35.7;  AE(90)a546 
eq.  (10)  with  a— 131,  b-582;  AE(90)>546 

eq.  (10)  with  a-305, 8-303;  AL(90)&  -335 

eq.  (10)  with  a-2.13  8-574;  AL(90)<  -335 

9 

eq.(8) 

«i(9) 

10 

eq.  (10)  with  a— 337,  b=503 

eq.  (10)  with  a-395,  8-23.7  . 

11 

eq.  (8a);  AE(90)S575 

eq. (11) withe— 103, d=.00134;  AE(90)>575 

eq.  (9a);  AL(90)&-450 

eq.  (11)  with e— 5.68, d— .00259;  AL(90)<-450 

12 

eq.  (8a);  AE(90)S547 

cq.  (11)  with  c— 11.7,  d«0.0013;  AE(90)>547 

eq.  (9a);  AL(90)i:-300 

eq.(ll)withc—m,d— .00193;  AL(90)<-300 

13 

eq.  (8a);  AE(90)S57S 

eq.  (11)  with  e— 530,  d-.00181;  AE(90)>575 

eq.(9a);  AL(90)a-450 

eq.(ll)withc—103,d— .00202;  AL(90)<-450 

14 

eq.(8) 

eq.  (9a);  AL(90)fc-450 

eq.  (11)  with  c-  -30.0,  d— .000983;  AL(90)<  -450 

15 

eq.  (10)  with  a=-317,  b=77.6 

eq.  (10)  with  a— .474, 8—283 

16 

eq.  (10)  with  a=-369,  b=40.4 

eq.  (10)  with  a— 351,  8-15.1 

17 

eq.  (8a);  AE(90)5547 

eq.  (10)  with  a= -.0585,8—28;  AE(90)>547 

eq.  (^);  AL(90)fc-344 

eq.  (10)  with  a-.080, 8-  -323;  AL(90)<  -334 

18 

eq.  (10)  with  a=-.102,  b=153 

eq.  (10)  with  a-.129,  8-7.75 

19 

eq.  (8) 

eq.  (9a);  AL(90)&-450 

eq.  (11)  with  c-  -333,  d-  -.000869;  AL(90)<  -450 

20 

eq.(8) 

eq.  (9a);  AL(90)&-450 

eq.  (11)  with  c— 29.1,  d— .00101;  AL(90)<-450 

21 

eq.(10)witha— .11,8—38;  AE(90)Sl205 

eq.  (10)  with  a— 1.9,8-2120;  AB(90)>1205 

eq.(10)witha-.lS5,b— 173;  AL(90)2:-1165 
eq.  (10)  with  a-301, 8-2242;  AL(90)<  -1165 

AL(90)  (nT) 


Fig.  29.  The  relationships  for  Dst  as  a  function  of  yLE(90)  (a)  andylL(90)  (b).  The 
general  relationships  from  equations  (8)  and  (9)  are  bold  and  the  relationships  for 
storms  1-21  are  thin. 


Table  5.  Correlation  coefficients  for  Dst  versus 
Dsr(AE(90))  and  £)sr(/4L(90))  for  storms  1-21. 
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Storm 

Carr.  Coelt  for 
Dsl(AE(90)) 

Carr.  Coell  for 
Dst(AL(90)) 

1 

0.711 

0.779 

2 

0J845 

0S95 

3 

0.711 

0.799 

4 

0.925 

0.960 

5 

0.897 

0S97 

6 

0.790 

0.905 

7 

0.745 

0.797 

8 

0.768 

0.884 

9 

0.931 

0.944 

10 

0.881 

0.879 

11 

0S83 

0.768 

12 

0.755 

0.799 

13 

0.919 

0S52 

14 

0.906 

0.897 

15 

0.795 

0S96 

16 

0.869 

0886 

17 

0.939 

0.943 

18 

0.903 

0888 

19 

0.892 

0.910 

20 

0.911 

0.931 

21 

0.777 

0816 
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CHAPTER  6 
DISCUSSION 

Results  of  Dst-AE  Comparisons 

The  results  of  Chapter  4  seem  to  show  that  during  storms  there  is  a  clear 
relationship  between  the  ring  current  and  electrojets.  Both  are  energized  by  the 
same  source  and  are  primarily  directly  driven.  While  most  researchers  use  the  1-hour 
AE  index  regularly  to  track  substorm  activity,  FFT  analysis  (see  Appendix  B)  shows 
(and  Fairfield  [1992]  agrees)  that  the  l-hour  AE  is  primarily  directly  driven,  especially 
since  the  actual  substorm  component  from  cross-tail  current  disruption  has  a  higher 
frequency  {Lee  et  aL  [1985];  Liu  et  al  [1988];  and  L.  Zhu,  private  communication, 
1993).  This  is  justified  by  the  result  of  AE^,^  (and  AL^)  giving  a  better  correlation 
than  AE  (and  AL)  since  the  smoothing  technique  was  shown  to  be  equivalent  to 
removing  higher  frequencies.  It  is  more  likely  that  the  1-hour  AE  index  contains 
primarily  the  directly  driven  component  of  the  auroral  currents  and  very  little  (if  any) 
loading-unloading  (substorm)  component  (this  view  is  supported  by  Fairfield  [1992]). 
This  would  explain  the  disagreement  of  those  who  claim  that  substorms  are  directly 
driven  {i.c.,Akasofu  [1981a])  by  comparison  with  the  solar  wind  (since  the  solar  wind 
components  used  drive  magnetospheric  convection)  and  those  who  claim  that  a 
substorm  has  both  directly  driven  and  loading-unloading  components  based  on 
magnetic  reconnection  and  plasmoid  observations  (i.e.,  McPherron  [1991]  and 
references  therein). 


59 


From  this,  then,  one  can  assume  that  when  the  IMF  turns  southward,  piasma 
flow  due  to  convection  increases  and  reaches  some  steady-state  value,  and  this  is 
reflected  in  the  elevated  ylE  level.  As  magnetic  reconnection  and  current  disruption 
occur  in  the  magnetotail,  substorms  occur  and  are  superposed  on  the  baseline  (as 
stated  by  Lee  et  al.  and  Du  et  al.)  but  are  of  little  additional  influence  on  the  1-hour 
AE  (see  Appendix  B  where  r  applied  to  the  1-minute  AE  is  roughty  equal  to  t 
applied  to  the  1-hour  XE),  so  the  1-hour  .4E  continues  to  reflect  the  directly  driven 
baseline  value.  At  the  same  time,  the  ring  current  is  being  energized  by  the  increased 
convection  and  is  also  direct^  driven.  This  continues  until  the  IMF  turns  northward 
and  the  enhanced  activity  returns  to  quiet  levels. 

With  this  picture  in  mind  the  results  of  using  the  t  weighting  on  the  AE  indices 
are  understandable,  since  l-hour  AE  and  Dst  are  both  measuring  a  directly  driven 
current  (driven  by  the  same  source)  that  differ  only  by  how  the  currents  respond  to 
injection  (this  is  also  supported  by  the  work  of  Kamide  and  Fukushima  [1971]).  The 
time-weighted  accumulation  method  is  equivalent  to  simulating  the  AE  current 
injection  as  an  injection  into  another  current  system  with  a  slower  response  time. 
The  result  is  that  the  new  indices  derived  from  t=0.9  bear  a  very  close  resemblance 
to  the  Dst  index,  indicating  that  the  theoretical  picture  outlined  above  has  some 
degree  of  authenticity. 

While  AC/(t)  did  not  correlate  well  to  Dst,  its  inclusion  in  AE(t)  did  not  cause 
a  significant  reduction  in  the  correlation  of  AE(t).  Although  it  could  be  argued  that 
AU  contributes  much  less  to  AE  and  is  not  as  much  of  a  factor,  it  could  also  be 
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argued  that  the  eastward  electrojet  closes  through  the  partial  ring  current  (Kamide 
and  Fukushima  [1971])  and  so  should  be  included  in  the  analysis  since  the  partial  ring 
current  contributes  to  Dst. 

Modeling  Dst  from  AE(t)  and  AL(t) 

Extending  the  above  comparisons  for  t=0.9  to  an  entire  year  of  data 
demonstrated  the  possibility  that  Dst  could  actually  be  modeled  or  forecasted  from 
AE(90)  and  AL(90).  The  fact  that  the  derived  Dst  followed  the  actual  Dst  very 
closely  at  low  values  shows  that  even  during  a  quieter  period  the  current  systems  are 
closely  tied  together  in  their  source,  indicating  that  the  above  relationship  between 
the  current  systems  is  a  permanent  one. 

Even  though  the  initial  comparisons  of  Dst  toAE(90)  andAL(90)  for  storms 
showed  some  similarity,  different  functional  relationships  between  Dst  and  AE(90) 
and  Dst  and  AL(90)  were  found  for  many  of  the  storms,  although  many  of  the 
relationships  were  very  similar  to  the  original  relationships  (8)  and  (9).  This  indicates 
the  possibility  that  some  other  physical  process  is  not  being  included  during  storm 
periods,  which  makes  the  relationships  different.  While  there  could  be  a  variety  of 
other  factors  to  consider,  an  initial  analysis  of  the  data  showed  no  dependence  on 
solar  cycle  or  season.  The  different  relationships  could  also  suggest  a  deficiency  in 
the  measurement  of  the  current  systems,  which  is  possible  given  all  the  inadequacies 
already  noted  in  Chapter  3. 
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Although  the  general  relationships  (8)  and  (9)  work  veiy  well  in  spite  of  the 
many  problems  with  the  indices,  it  is  useful  to  speculate  on  the  reasons  why  they  do 
not  work  better: 

1)  While  the  partial  ring  current  may  or  may  not  be  reflected  in  xhtAU  index, 
it  no  doubt  has  a  different  time  constant  (as  shown  by  Kamide  and  Fukushima)  than 
the  symmetric  ring  current.  It  would  therefore  reduce  the  relationship  between  Dst 
and  AE(f)IAL(T)  using  a  r  that  fits  the  symmetric  ring  current. 

2)  During  the  first  few  hours  of  ring  current  energization,  the  ring  current 
cannot  be  symmetric;  therefore,  that  assumption  about  the  Dst  index  is  no  longer 
valid. 

3)  The  AE  stations  are  not  always  in  position  to  measure  the  fuU  intensity  of 
the  electrojets  so  that  measurements  are  missing  that  would  change  the  character  of 
AE{t). 

4)  The  application  of  t  to  the  AE  indices  means  that  one  time  constant  is 
chosen  and  remains  constant,  while  most  ring  current  researchers  have  established 
that  the  time  constant  of  ring  current  decay  changes  during  the  course  of  most 
storms. 

5)  The  effect  of  magnetospheric  compression  is  not  removed  from  Dst  so  that 
the  ring  current  is  not  always  the  only  current  being  measured. 

6)  It  is  possible  that  electrojet  effects  are  included  in  the  ring  current, 
especially  during  intense  storms  when  the  electrojet  moves  equatorward. 
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7)  It  is  possible,  as  Akasofu  et  at  [1983]  stated,  that  AE  becomes  saturated 
during  intense  storms  and  is  not  really  representing  the  electrojets,  although  no 


evidence  of  this  is  found. 


CHAPTER  7 


CONCLUSIONS  AND  FURTHER  RESEARCH 

Conclusions 

The  conclusions  that  can  be  drawn  from  the  discussion  in  the  previous  chapter 
are  summarized  as  follows: 

1)  The  1-hour^  indices  are  totally  (or  nearly  so)  directly  driven  in  that  they 
contain  almost  no  contribution  from  the  substorm  current  wedge.  This  would  resolve 
the  apparent  conflict  between  those  researchers  who  argue  that  substorms  are  directly 
driven  versus  those  who  say  that  substorms  have  both  a  directly  driven  and  a  loading¬ 
unloading  component. 

2)  A  time-weighted  accumulations  method  applied  to  the  1-hour  AE  indices 
gives  a  high  correlation  to  1-hour  Dst  during  intense  storms. 

3)  The  time-weighted  accumulation  results  show  a  close  relationship  between 
the  ring  current  and  the  electrojets,  with  a  difference  of  only  response  time  to 
injection. 

4)  Most  of  the  time  there  is  a  clear  one-to-one  relationship  between  Dst  and 
X£(90)/4L(90).  During  storms  this  relationship  may  change  but  one  still  exists.  The 
reason  for  the  changing  relationship  remains  to  be  explored. 

5)  The  possibility  exists  of  modeling  Dst  in  real  time  iiAE  data  is  available. 
The  U.  S.  Air  Force  is  currently  developing  a  chain  of  magnetometer  stations  which 
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would  provide  real-time  IQ)  data  with  15-minute  resolution.  This  could  easily  be 
configured  to  provide  an  estimate  of  theXE  indices  from  the  magnetometer  data. 

The  main  thrust  of  magnetospheric  research  for  the  past  25  years  has  been  in 
the  area  of  substorms,  with  the  belief  that  this  would  in  turn  lead  to  a  better 
understanding  of  storms  (i.e.,  the  ring  current).  These  results  seem  to  suggest  that 
this  notion  may  be  somewhat  misguided  in  that  substorms  make  very  little 
contribution  to  the  overall  auroral  current  ^tem,  which  has  already  been  greatly 
enhanced.  These  enhanced  background  currents  are  the  ones  that  show  a  direct 
relationship  to  Dst  (and  hence  the  ring  current)  during  a  storm,  while  the  short-term 
fluctuations  from  the  substorm  current  wedge  are  left  out  of  the  picture  entirely. 

Areas  For  Further  Research 

Some  possible  areas  to  pursue  in  connection  with  this  research  are: 

1)  Finding  the  physical  mechanism  responsible  for  the  different  functional 
relationships  between  Dst  and  AE(90)/AL(90)  during  intense  storms. 

2)  Adjusting  t  to  maximize  correlations.  Although  r=0.9  was  found  to  give 
the  best  overall  results,  different  values  of  t  worked  better  for  some  storms.  A 
different  t  for  each  storm  would  undoubtedty  produce  better  correlations. 

3)  Varying  t  throughout  the  storm  to  maximize  correlations.  This  would  be 
an  extension  of  the  previous  topic  and  would  parallel  the  work  of  most  researchers 
who  model  Dst  fi'om  IMF  data  and  use  a  variable  decay  parameter  (i.e.,  the  methods 
reviewed  by  Feldstein  [1992]). 
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7)  It  is  possible,  as  Akasofu  et  at  [1983]  stated,  that  AE  becomes  saturated 
during  intense  storms  and  is  not  really  representing  the  electrojets,  although  no 
evidence  of  this  is  found. 
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CHAPTER? 

CONCLUSIONS  AND  FURTHER  RESEARCH 

Conclusions 

The  conclusions  that  can  be  drawn  from  the  discussion  in  the  previous  chapter 
are  summarized  as  follows: 

1)  The  1-hour^  indices  are  totally  (or  nearly  so)  directly  driven  in  that  they 
contain  almost  no  contribution  from  the  substorm  current  wedge.  This  would  resolve 
the  apparent  conflict  between  those  researchers  who  argue  that  substorms  are  directly 
driven  versus  those  who  say  that  substorms  have  both  a  directly  driven  and  a  loading- 
unloading  component 

2)  A  time-weighted  accumulations  method  applied  to  the  1-hour  AE  indices 
gives  a  high  correlation  to  1-hour  Dst  during  intense  storms. 

3)  The  time-weighted  accumulation  results  show  a  close  relationship  between 
the  ring  current  and  the  electrojets,  with  a  difference  of  only  response  time  to 
injection. 

4)  Most  of  the  time  there  is  a  clear  one-to-one  relationship  between  Dst  and 
.(4£(90)44L(90).  During  storms  this  relationship  may  change  but  one  still  exists.  The 
reason  for  the  changing  relationship  remains  to  be  explored. 

5)  The  possibility  masts  of  modeling  Dst  in  real  time  if  AE  data  is  available. 
The  U.  S.  Air  Force  is  currently  developing  a  chain  of  magnetometer  stations  which 
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APPENDIX  A.  STORM  PLOTS 
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Fig.  A-1.  Plots  for  storm  1.  Dst  (bold)  is  plotted  versus  AE,AE^AE(t=0.75),  and 
AE(t=0.9)  (a).  The  same  is  done  for  AL  (b),  AU  (c),  and  AO  (d). 
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Fig.  A-2.  Plots  for  storm  2.  Dst  (bold)  is  plotted  versus  AE,AE^.4£(t=0.75),  and 
A£(t=0.9)  (a).  The  same  is  done  for  AL  (b),  AU  (c),  and  AO  (d). 
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Fig.  A-3.  Plots  for  storm  3.  Dst  (bold)  is  plotted  versusvlE,  AE(t=0.75),  and 
AE(t=0.9)  (a).  The  same  is  done  for  AL  (b),  Af/  (c),  and  AO  (d). 
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Fig.  A.5.  Plots  for  storm  5.  Dst  (bold)  is  plotted  versus  AE,y4£^AE(T=0.75),  and 
AE(t=Q.9)  (a).  The  same  is  done  ioiAL  (}3\AU  (c),  and  ^<7  (d). 
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Fig.  A.8.  Plots  for  storm  8.  Dst  (bold)  is  plotted  versus  AE,XE^AE(t=0.75),  and 
AE(t=0.9)  (a).  The  same  is  done  for  AL  (b),  Af/  (c),  and  AO  (d). 
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Fig.  A9.  Plots  for  storm  9.  Dst  (bold)  is  plotted  versus  AE,./4£^.(4£(t=0.75),  and 
XE(t=0.9)  (a).  The  same  is  done  for  AL  (h),AU  (c),  and  AO  (d). 
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Fig.  A.12.  Plots  foi  storm  12.  Dst  (bold)  is  plotted  versus  A£,XE^AE(t=0.75),  and 
XE(t=0.9)  (a).  The  same  is  done  ioiAL  (h),AU  (c),  andylO 


Fig.  A.13.  Plots  for  storm  13.  Dst  (bold)  is  plotted  versus  >4£(t=0.75),  and 

AE(t=0.9)  (a).  The  same  is  done  for  AL  (b),  AU  (c),  and  AO 
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Storm  15 

(a)  (b) 


Fig.  A.15.  Plots  for  storm  15.  Dst  (bold)  is  plotted  versus  XE,y4£^AE(T=0.75),  and 
AE(t=0.9)  (a).  The  same  is  done  for  AL  (b),  AU  (c),  and  AO 
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Fig.  A- 16.  Plots  for  storm  16.  Dst  (bold)  is  plotted  versus 
AE(t=0.9)  (a).  The  same  is  done  for  AL  (h),AU  (c),  and  AO  ^ 
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Fig.  A.17.  Plots  for  storm  17.  Dst  (bold)  is  plotted  \tTsmAE,AE^AE(T=Q.15),  and 
>4£(t=0.9)  (a).  The  same  is  done  for  AL  (b),  AU  (c),  and  AO  (5. 
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Fig.  A.20.  Plots  for  storm  20.  Dst  (bold)  is  plotted  ytT%\3&AE,AEg^AE{T=Q.15),  and 
AE(t=0.9)  (a).  The  same  is  done  ioxAL  (p\AU  (c),  snd  AO  (d). 
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Fig.  A.21.  Plots  for  storm  21.  Dst  (bold)  is  plotted  versus  ylE,AE^XE(T=0.75),  and 
AE(t=0.9)  (a).  The  same  is  done  for  AL  (b),AU  (c),  and  AO  (d). 
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FOURER  ANALYSIS  OF  1  MINUTE  AE  DATA 


A  Fourier  Transform  analysis  was  done  on  one  sample  storm  to  determine  the 
properties  of  the  1-minute  AE  index  as  it  relates  to  the  1-hour  index  and  the  various 
manipulations  performed  on  it.  The  storm  selected  was  storm  11,  July  13-16,  1982 
because  it  is  fairly  representative  of  a  normal  storm. 

First,  the  Fourier  coefficients  were  obtained  from  the  following  relation 
[Chatfield,  1975]: 


a^=  X 


(Al) 


where  N  is  the  number  of  datapoints  and  p  ranges  from  1  to  N/2-1.  The  1-minute 
index  was  then  reconstructed  using: 


2 

^  [a^cos(^^)  +  f»^sin(^^)]  +  a^cos(nO 


(A2) 


with  the  higher  frequencies  removed.  Figure  A.22  shows  that  when  frequencies 
greater  than  69  MHz  are  removed  (corresponding  to  periods  of  4  hours  and  lower). 
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the  result  is  very  similar  to  the  1-hour  AE  index.  When  frequencies  higher  than  28 
MHz  are  removed,  it  is  very  similar  to  achieved  by  space-centered  averaging 
of  the  1-hour  AE  index  with  10  iterations.  This  demonstrates  that  the  1-hour  AE 
index  contains  only  lower  frequency  components  of  the  electrojet  variations  and  that 
AEg^^  contains  only  increasingly  lower  components  as  the  number  of  iterations  is 
increased. 

Next,  time-weighted  accumulations  were  applied  to  the  1-minute  AE  index. 
In  Figure  A.23  it  can  be  seen  that  a  t  of  0.9883  (corresponding  to  a  decay  time  of 
-10  hours  for  1-minute  data)  applied  to  the  1-minute  AE  is  virtually  identical  to  a  t 
of  0.9  (corresponding  to  a  decay  time  of  -10  hours  for  1-hour  data)  used  on  the  1- 
hour  AE.  Similar  results  are  found  when  t=0.9883  is  used  for  1-minute  AE  with 
higher  frequencies  removed.  Therefore,  when  time-weighted  accumulations  are  used, 
very  little  is  lost  when  using  the  1-hour  AE  index  as  opposed  to  the  1-minute  index. 


Fig.  A.22.  One  minute  yl£  (top)  and  1  minute  AE  with  higher  frequencies  removed 
compared  to  1  hour  AE  (middle)  and  (bottom). 


Fig.  A.23,  Time-weighted  accumulations  of  1  minute  yl£  (t=0.9883),  both  unaltered 
and  with  higher  frequencies  removed,  compared  to  time-weighted  accumulations  of 
1  hour^^E  (t=0.9). 
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DST  MODELING  FROM  SOLAR  WIND  PARAMETERS 


As  stated  by  Feldstein  [1992]  in  his  review  of  ring  current  modeling,  the 
method  outlined  by  Pisarsky  et  at  [1989]  is  one  of  the  best  for  deriving  Dst  from  IMF 
data,  in  order  to  give  examples  of  the  current  state  of  Dst  modeling,  this  method  will 
be  applied  to  those  storms  of  the  21  that  have  as  complete  a  set  of  IMF  data  as 
possible. 

Pisarslgr  et  al.  used  this  relation: 


DR^(2--)  +  2F^ 


2  + 


(A3) 


where  DR  is  the  measure  of  the  ring  current  (Dst  minus  magnetospheric 
compression)  and  F  is  the  ring  current  injection  function  modeled  by: 

F  =  (8,2xl0-^)v(B^-0.67a)-14.1xl0-’(v-300)+9.4, 
for  v(B^-0.67o)  <  -1146  (injection) 

F  =  (-14.1  xlO-’(v- 300),  (A5) 

for  v(B-0.61a)  >  -1146  (decay) 


where  a  is  the  IMF  modulus.  The  decay  parameter  is  obtained  from: 


=  2.4  +  13«®°*^  (main  phase) 

=  10.0  +  1.84  (recovery  phase) 


(A6) 
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using  the  F  or  DR  value  for  that  hour.  Since  no  definition  of  the  IMF  modulus  a 
could  be  found,  it  was  set  to  0.  The  results  of  applying  this  method  to  storms  2,  3, 
4,  6,  and  9  are  shown  in  Figures  A.24-A.28. 


Storm  2 


Fig.  A-24.  Dst  (bold  line)  compared  to  Dst  derived  from  the  method  of  Pisarsky  et  al. 
for  Storm  2. 


Fig.  A.25.  Dst  (bold  line)  compared  to  Dst  derived  from  the  method  of  Pisarslg^  et  al. 
for  Storm  3. 


Storm  4 


0  so  100 


Hours 

Fig.  A.26.  Dst  (bold  line)  compared  to  Dst  derived  from  the  method  of  Pisarslqr  et  al. 
for  Storm  4. 


102 


Storm  6 


Fig.  A.27.  Dst  (bold  line)  compared  to  Dst  derived  from  the  method  of  Pisarsky  et  al. 
for  Storm  6. 


Storm  9 


Fig.  A.28.  Dst  (bold  line)  compared  to  Dst  derived  from  the  method  of  Pisarsky  et  al, 
for  Storm  9. 
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